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Abstract: We developed a predictive model describing harmonic gener- 
ation and intermodulation distortions in semiconductor optical amplifiers 
(SOAs). This model takes into account the variations of the saturation 
parameters along the propagation axis inside the SOA, and uses a rigorous 
expression of the gain oscillations harmonics. We derived the spurious-free 
dynamic range (SFDR) of a slow light delay line based on coherent 
population oscillation (CPO) effects, in a frequency range covering radar 
appUcations (from 40 kHz up to 30 GHz), and for a large range of injected 
currents. The influence of the high order distortions in the input microwave 
spectrum is discussed, and in particular, an interpretation of the SFDR 
improvement of a Mach-Zehnder modulator by CPOs effects in a SOA is 
given. 
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1. Introduction 

A continuously tunable optical delay line is a key element for a large variety of microwave 
photonics applications, including the control of optically fed phased array antennas, the filter- 
ing of microwave signals or the synchronization of optoelectronic oscillators |[Tl|2l[3|. Slow 
light based tunable delay lines have been therefore intensively studied over the past few years. 
The main focus has been to understand and modelize the underlying phenomena H, in order to 
increase the delays and/or the bandwidth of the studied effect Among the different slow 

and fast light architectures, coherent population oscillations (CPOs) in a Semiconductor Optical 
Amplifier (SOA) offer attractive operational advantages in terms of compactness, integrability 
and possible parallelism, as well as a continuous tunability of the delay, or phase shift, through 
the injected current. Moreover, a recent demonstration of 360° phase shift obtained at up to 19 
GHz using this technology, makes it very promising towards integration in radar systems. 
The principle of slow light in SOAs can be understood as follows. An optical carrier with a 
sinusoidal modulation envelope that propagates in a SOA induces a carrier population oscilla- 
tion, via gain saturation in the semiconductor material. Since gain oscillations are in antiphase 
with the modulation envelope, and basically occur for modulation frequencies below the inverse 
of the carrier lifetime, they induce a dip in the RF gain spectrum of the SOA. Consequently, 
according to the Kramers -Kronig relations, the gain dip is associated with a large and positive 



refractive index dispersion, and hence a large group index, thus slowing down the sinusoidal 
modulation propagation velocity. This property, as well as corresponding phase and amplitude 
change of the modulation envelope, have been extensively studied, both experimentally or theo- 
retically |[8ij9][T0][TT]. In particular, this has provided the RF transfer function of the SOA-based 
delay line, similarly to any two-port electronic device. 

However, within the scope of an integration in a realistic radar system, it is also required to study 
the impact of these slow light architectures on the whole microwave photonics link, and in par- 
ticular the generation of spurious signals by harmonic and intermodulation products. One of the 
key parameters of a microwave photonics link for radar applications is indeed the spurious-free 
dynamic range (SFDR) |T2|. It represents the operational dynamic range of the link, i.e., the 
range of input (or output) RF powers, at a given frequency, for which all spurious signals at har- 
monic and intermodulation frequencies are below the noise level at the output of the link. The 
SFDR requirements for radar microwave photonics links are more or less stringent depending 
on the application, but are basically in the range of 105dB/Hz^/^ (ground based antennas) to 
1 15dB/Hz^/^ (airborne antennas). Most of the previously reported work on nonlinear frequency 
mixing in SOAs has been carried out in the frame of Radio-over-Fiber applications, or SOA- 
based in-line photodetection ifTSl [T4l [TSl . The experimental and operational context of these 
works significantly differs from radar applications. As an example, in |T3], the SOA is used 
for multifrequency signal up-conversion, where incoherent frequency channels are considered. 
However, a typical radar microwave photonics link involves a single electro-optic modulator 
that indistinctly transfers all the spectral components of the received RF signal (including a 
jammer, for instance) to the optical carrier Coherent beating between spectral components can 
thus not be avoided, and the SFDR measurements reported in IfTSl can therefore not be extrap- 
olated to slow light based tunable delay lines for radar systems. 

In this paper, we first present a predictive model to determine the harmonic generation level, 
which takes into account the spatial evolution of the saturation parameters along the SOA 
length. We also discuss the expression of the gain modulation harmonics that is used in com- 
mon models, and show that it can lead, under certain conditions, to a significant error on the 
harmonic's estimate. Then, we generalize our model to end up with the third order intermodula- 
tion distortion (IMD3), and compute the IMD3 level and the SFDR over the frequency range of 
interest for radar applications (up to 20GHz), and for a large range of injected currents. Finally, 
we discuss the influence of the input optical spectrum on the dynamic range of the delay line by 
comparing the simulation results obtained for an ideal spectrum, i.e., a perfect sinusoidal modu- 
lation of the optical beam at the input of the SOA, and for a realistic spectrum, namely, when the 
RF signals to be delayed are transferred on the optical carrier using a standard Mach-Zehnder 
modulator 

2. Principle of harmonics calculation 

We consider an optical field E{z,t) which propagates along a traveling wave semiconductor 
optical amplifier (SOA). The interaction of light with carriers in the SOA is governed by the 
well known carrier rate equation and field propagation equation HI: 



dN{z,t) 
dt 
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where N{z,t) is the carrier density, / the current injected in the SOA, Ty the carrier lifetime, q 
the elementary electric charge, V the SOA active volume, g the optical gain, F the confinement 



factor, and y\ stands for the internal losses. In order to find the level of the generated harmonics, 
we first consider that the input optical field is modulated at the RF frequency H. § and 
are hence all time-periodic functions with a fundamental frequency of i2. They can therefore 
be written into Fourier harmonic decompositions: 
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where N{z) and g{z) respectively denote the DC components of the carrier density and of the 
optical gain. a{z) is the SOA differential gain, defined as a{z) ~ dg/dN. Defining gi^ as the 
oscillating component of the gain at frequency kO., and considering only a finite number K of 
harmonics, the carrier rate equation ([T]) can be written in an equivalent matrix formulation: 
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where = /v(l +Mo//.v — ikD.Ts), and ao = Mq is the DC optical intensity. denotes the 
saturation intensity and is defined as /j = ho/aZs- It is worth mentioning that is obtained at 
the first order of equation ([T]), when mixing terms are not considered. One can also notice that 
gl = g_i^. For the sake of clarity, we did not write the dependence on the propagation coordinate 
z in (|6]l. It is however important to note that in the following, .^V, g, a, Zs, Is, and consequently the 
a^'s are all actually functions of z- Their variations along the propagation axis is then taken into 
account, unlike most of the reported models in which effective parameters are used lISl fTOlfTTI . 



2.1. Small RF signal formulation 

In order to solve the propagation equation one first has to determine the expressions of 
the gain harmonic components g^-'s, as functions of the harmonic components of the optical 
intensity M^^'s, of g, and of /y and Tj. Under small RF signal approximation, i.e. considering 
a small modulation index of the optical carrier, one can assume that |Mj^_i | 3> |Mi.| for A: > 
(and |Mi._i | <C |M<.| for k < 0). It can also be noticed that the same relations hold for the N^'s 
since the carrier density oscillations are induced by the illumination oscillations through gain 
saturation. On the basis of this assumption, the commonly used expression for the coefficients 
gkisMMM: 



where 



g = TicoiI/qV-N/T,)/Mo. (8) 

In other words, this is equivalent to consider as nonzero only the diagonal and center column 
of the matrix in However, the quantities to be compared in the inversion of (|6]l are not the 
|M^|'s, but the terms |Mj^_p x gp\- Consequently, the hypothesis {Mi^^i \ ^ \Mit \ leads to neglect 
only the terms for which \p\ + \k — p\ > \k\, the others being of the same order Under these 
conditions, ^ is reduced into; 
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The resolution of this equation gives identical expressions for g and g±\ as in (|7]) and ([8]), 
indicating no change on the fundamental component of the output modulation, compared to 
common models. However, for |^| > 1, additional terms appear, whose importance will be 
emphasized later on. The expressions of g2 and ^3 are hence: 
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Once again, /, and t, in equations (|7]) to ( fTTT ) are functions of the propagation coordinate z, 
through their dependence on the static carrier density N{z). Doing so, the relations between /,, 
Tv and the physical device constants are preserved, which ensures the predicting capability of 
the model when the operating conditions (optical input power or bias current) are changed. As 
reportes in |fT9l , a simple measurement of the unsaturated gain as a function of the SOA bias 
current, gives a first relation between g and N /is'- 



g = a- 
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(12) 



where a and j3 are empiric coefficients determined experimentally by measuring the SOA's 
small signal gain. For our simulations, we used a = 5.88 • lO^m^' and )3 = —1.84 • lO^^m^.s^', 
corresponding to the COVEGA InGaAsP/InP quantum well SOA available in our laboratory 
IIT9I . Then, solving the system constituted of (|8]l and (fT2] i. we obtain g{z) and N/ T,(z) as func- 
tions of a, P, I and the local DC optical intensity Mo(z). Finally, we model the carrier lifetime 



in our SOA using the well known expression: 
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where A, B and C are the carrier recombination coefficients of the semiconductor structure, 
and correspond respectively to the non-radiative, spontaneous, and Auger recombination coef- 
ficients. Equation ( fTSl ) and the expressions of g{z) and N /Ts{z) then enable to derive T, and /, 
as functions of /, Mq{z) and A, B and C. 

According to lfT9l . A, B and C are the only adjustment parameters of the model: they are de- 
termined by adjusting the simulated and measured fundamental RF transfer functions of the 
SOA. In the case of the COVEGA SOA we consider in this paper, we obtainedA = 2- lO'^s-i, 
B = 1.2 ■ 10^"^cm^.s^\ C = 1.8 • lO^-^'cm^.s^'. It is worthwhile to mention that these values 
were found to be valid for the full range of injected currents or optical input powers, proving 
the predictive capability of this approach. 



2.2. Large RF signal formulation 

In the case of a large modulation index, the previous assumption |M^^_i | jM^^j do no longer 
apply. Equation ^ has to be rigorously solved. Moreover, equation (O is no longer valid, 
and consequently, g. Is and Tj cannot be obtained as in the small signal case. We thus use the 
following iterative procedure: in a first step, we substitute N/Ts, It and in (|6]l by their small 
signal values N/ Ti**', and ri"' as obtained in section ITT] The central matrix of (|6]l, referred 
as D, can then be inversed. The coefficients of D 
then becomes: 



are denoted (5^,^). Equations (|7]i and dH) 
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Similarly to the small signal case, using equations ( fT2b . ( fTsT i and ( fTSl ). we obtain .^V/ri^', /J'' 
and Ty'' as functions of /, A, B, C and Mi^{z). This procedure is repeated until convergence 
of N/Ts"\ /i"' and ri"', which typically occurs after a few tens of iterations. The propagation 
equation ^ can now be expressed in a matrix formulation similarly to ifTTI and numerically 
solved. 
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2.3. Comparison with reported models 

To demonstrate the importance of the additional terms in the right hand side of the expressions 
of g2 and ^3 in equations (fTol i and (fTTT i. we calculated the third harmonic power i/3 at the output 



of a 1 .5mm-long SOA. The optical carrier is modulated using a Mach-Zehnder modulator The 
initial conditions for the resolution of equation ( fTSI l are hence 1201 : 

Mo,,„=P,„ 
M\.in = -Pi„xJi{m) 
< M2,,„ = , (17) 

where m is the modulation index, and 7^ denotes the A:''' order Bessel function of the first kind. 
For this example, we chose a modulation index of 0.01, which ensures the validity of the small 
RF signal hypothesis. 



SOA parameter 


Value 


Facet to facet linear gain 


37.6dB @I=500mA 


Saturation power 


15dBm @I=500mA 


Linear losses 


13.6cm-' 


Coupling losses 


3.3dB/interface 


Amplifier length 


1.5mm 



Table 1. List of studied SOA parameters 



On the other hand, we measured experimentally the third harmonic level at the input and at 
the output of a commercial SOA, whose main physical parameters are listed in Table [T] These 
parameters were also used in the simulations. It is worth mentioning that the saturation power 
given in the table is the measured output power for which the unsaturated optical gain is reduced 
by 3dB. Therefore it does not correspond to Is in the model, which is a local parameter, varying 
along the propagation axis. The optical input power is lOdBm, for which the optical gain is 
strongly saturated, and which leads to the maximum RF phase shift that we obtained with this 
device [19] . Equation ( fTSI l is solved using the expression of the g^'s given first by the common 
equations (|7|i and ([S]), and then by equations ( fTol i and (fTTT i containing the additional terms. 
Figure |l(a)| shows both experimental (red circles) and theoretical (solid lines) third harmonic 
photodetected power, normalized to the third harmonic power at the input of the SOA. The 
theoretical curves are evaluated according to: 

//3-27Jt72^|M3,„„,x5|2 (18) 

where R and rjpi, are respectively the photodiode resistive load (50i2) and efficiency (equal to 
0.8). S denotes the SOA modal area. 

As can be seen on figure[T(a)j for a modulation frequency Q. above AGHz, the two models are 
equivalent, and the third harmonic power naturally tends to the value it would have if no CPO 
effects were present (i.e. M3 „,„ is equal to M3 ,„ multiplied by the optical gain Gopt)- However, 
at low frequencies, the difference between the two models reaches 20dB. The experimental 
measurement shows a very good agreement with our model including the additional terms, 
which confirms the vaUdity of our approach. The large discrepancy with common models can 
be qualitatively understood considering the asymptotic case where H tends to zero. Using equa- 
tions ( [Tot and (fTTT i in the propagation equation ( fTSI l. one obtains for the third harmonic term. 




Fig. 1. (a) Third harmonic generation as a function of the RF frequency. Green solid hne: 
common expressions of gk (equations Q and ^) are used. Red solid line: use of the 
rigorous expressions of gk (equations dlOb and The red circles represent experimental 
measurements, (b) Asymptotic case when f2 — > 0: Evolution of terms A and B along the 
SOA, and third harmonic level calculated according to ( |19b (red solid line) and i20\ (green 
solid line). 



with Ht, < 1: 
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whereas, using the standard models based on equation O, one has: 
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Basically, the right-hand side of these two equations consists in three terms. The first two re- 
spectively describe the linear losses and amplification of M3. The third one, respectively de- 
noted A and B for equations (T% and ( l20b describe the energy transfer from Mi and M2 to M3 
through CPO process. Both A and B have a negative sign, which is consistent with the CPO's 
oscillations being in antiphase with the illumination oscillations. In figure [T(b)| are plotted the 
evolution along the propagation axis inside the SOA of terms A and B, and of M3, calculated ac- 
cording to equation iT% in red, and according to equation ( |20l ) in green. These results indicate 
that \B\ is much larger than \A\, and hence that conventional models overestimates the third har- 
monic generation term due to CPO for the low frequency part of the spectrum. Consequently, 
CPO contribution largely dominates the amplification term, and leads to a relatively high output 
level of third harmonic. On the other hand, when the gain harmonics are rigorously derived, one 
ends up with the term A, with a smaller magnitude which balances the amplification process. 
The resulting output Hj, level is therefore significantly reduced, and can even be under the input 
level, if the optical gain is saturated enough such as the overall amplification/generation terms 
are below the linear losses attenuation. 

This shows the importance of considering accurate expressions for gi^ even in the small signal 
situation, and especially for the low frequency part of the microwave spectrum. When Q. in- 
creases, since the 2'"' and 3"' terms of ( fTTI ) respectively evolve as Q.^ and Q.^, then equation 



(fTTl i progressively tends to the common expression given by equation This can be seen on 
figure |l(a)[ where the two curves finally coincide for Q. > AGHz- Consequently, the present 
analysis does not question the results obtained for instance in ITSll . where harmonics calcula- 
tions are carried out at a fixed frequency of IQGHz- However, as confirmed by experiments, 
a rigorous calculation is mandatory when the modulation frequency lies in the spectral region 
1—3 GHz where most of ground radars operate. 

3. Intermodulation distortion 

Intermodulation distortion (IMD) calculation is slightly different from what has been discussed 
in the above section. Indeed, the number of mixing terms that must be taken into account is sig- 
nificantly higher For radar applications a typical situation where the IMD plays a crucial role is 
that of a radar emitting at a RF frequency i2 1 , and facing a j ammer emitting at , close to il i . 
Both £l\ and ^2 are collected by the antenna and transferred to the optical carrier through a sin- 
gle electro-optic modulator The point is then to determine the nonlinear frequency mixing due 
to the coherent population oscillations (CPO) inside the SOA. In particular, the mixing prod- 
ucts at frequencies £^2 — £l\ (or £2i — £^2) and 2Q.2 — £2i (or 2£l\ — Q.2) — respectively called 
second (IMD2) and third (IMD3) order intermodulation distortions — have to be evaluated at 
the output of the SOA. 

The main difference with harmonic calculation is that the optical intensity, and hence the SOA 
carrier density A^, and the SOA gain g are no longer time-periodic functions of period £2, but of 
period 5£l = £22 — iij. Then, using a Fourier decomposition of jEp, A^, and g, one has: 

|£(z,f)P = £ M,(zK''^^"', (21) 
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We consider a typical radar frequency £2i of lOGHz, and a frequency spacing 5£2 of lOMHz. 
Here, for intermodulation distortion calculation, we assume that only the spectral components 
at £2i,2, 2£2i,2, and all their first order mixing products significantly contribute to the generation 
of IMD2 and IMD3, as illustrated in figure |2] The A^'s and the g^^'s are then reduced in 19 
elements vectors, and gathered into blocks, denoted Mhiock.j and ghiock.j- The block contains 
the mixing products with frequencies close to j xQ.^. According to the notation of figure |2] the 
carrier rate equation ^ and the propagation equation (|2]i can respectively be written as: 
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where D, y and Hi j are sub-blocks of the complete matrices of equations (|6|l and (fT6] l respec- 
tively. They describe the mixing of the spectral components around the and the /'' harmon- 
ics. From now on, the full procedure described in section|2]can be applied in the same iterative 
way to determine the ^^^'s. Is and T,, and to finally numerically solve the equation dZSl ). Similarly 
to equation (fTsl l. the photodetected RF power at 2^2 — Hi is then calculated through: 
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Fig. 2. Set of significant spectral components of A' and g, and associated index k in 
their Fourier decompositions, n is defined such as flj = n5Q.. 



4. Simulation results 



In the general situation depicted in section |3] the optical intensity at the output of a chirp-free 
Mach-Zehnder modulator, that is, at the input of the SOA, will be of the following form | 



^opt.w 



1 + cos(m[cos{D.it) + cos{D.2t)) + 



(27) 



where m is the modulation index, and (j) is the modulator phase bias. When developed into 
first kind Bessel functions, and according to the notation of figure |2] the corresponding input 
distribution of the M^^'s, at quadrature bias = 7r/2), is: 
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In particular, one can notice that the optical intensity at the input of the SOA contains a term at 
2^1 2 — ii2,i, denoted IMD3", and equal to J\ {m)J2{m). In a general way, when considering the 



propagation equation in its matricial form (eqs. ( fT6b and dZSTl), the propagation of the Fourier 
compounds of an optically carried microwave signal into the SOA can be seen as resulting 
from an amplification process (the diagonal terms of the matrix) and a generation process by 
frequency mixing through coherent population oscillations (non-diagonal terms of the matrix). 
In order to better understand the combination of these two effects in the case of IMD3 prop- 
agation, we conducted our simulations considering two initial conditions, with and without 
IMD3". The latest case corresponds to the use of a perfectly linear modulator, or as in 1 13J for 
Radio-over-Fiber applications, where modulations at iii and £22 are produced by two distinct 
modulators on two incoherent optical carriers. An alternative way to avoid IMD^f would be to 
use a Single-Side-Band modulator, as proposed in 112111221 . In a tunable delay line based on 
slow-light effects in a SOA, the optical group delay can be typically tuned either by varying the 
optical input power, or the SOA bias current. Although both situations can be easily simulated, 
in this paper, we only present the latter one, using the bias current, which is the most suitable 
for an implementation in a real radar system. Figure [3(a)| represents the simulated IMD3, cal- 
culated according to equation ( |26] |. as a function of the modulation frequency £2i, and for a set 
of SOA injected currents, ranging from 75 mA, slightly above the transparency (50 mA), up to 
the maximal current of 550mA. IMD3 is known to scale with {P'l^f)^, where Pj^p the RF power 
applied to the Mach-Zehnder modulator. Moreover, P^^p is proportional to m^. Thus, in order 
to get rid of the dependence on the modulation index, we chose to normalize IMD3 by m^. In 
figure [3(a)| the simulation results are plotted for m = 0.01 and m = 0.5. As expected, the curves 
are perfectly coincident. 




Fig. 3. (a) Computed third order intermodulation power IMD3, normalized to m^, as a func- 
tion of the modulation frequency, for various injected currents, (b) Corresponding phase of 
the beat-note at 2Q.2 — ■ Dashed lines: case of a perfectly linear modulator (IMD™ = 0). 
Solid lines: actual Mach-Zehnder modulator (IMD™ ^ 0). 

The expected low-pass filter behavior of the SOA is observed on figure [3(a)] when no beat- 
note at 2Q.2 — £2i is injected into the SOA (dashed lines) lfT6l . The situation is completely 
different when the nonlinearities of the modulator used to transfer the RF signal to the optical 
carrier are taken into account in the simulations (solid lines). First, as for the third harmonic 
generation (see section |273] l, the asymptotic values for the high frequency part of the spectrum 
naturally corresponds to IMDj" multiplied by G^p,, for the same reason as the low-pass be- 
havior mentioned above, i.e. that CPO effects are roughly restricted to frequencies below 1 /Tj. 
Conversely, for the low frequency part of the spectrum, the intermodulation products can be 
seen as the result of the combination of two effects, namely the amplification of the incident 
beat-note at 2Q.2 — £2i, and the generation of a 2Q.2 — £2i mixing term due to the gain mod- 



ulation induced by the CPO effects inside the SOA. The major difference between these two 
effects is that the amplified beat-note is roughly in phase with the incident one, whereas for the 
CPO effects, due to the tt -phase shift between the intensity and gain modulations, the phase of 
the CPO-induced beat-note at 2Q.2 — iii is also ;r-shifted with respect to the incident beat-note. 
To better understand the results of figure |3(a)| and particularly the dip observed in the IMD3 
curves, we also computed the evolution of the phase of the beat-note term M2n2-ni versus mod- 
ulation frequency (see figure [3(b)] i. At low frequencies below I/Tj, CPO effects dominates, as 
indicated by the phase curves distribution around —n, whereas at high frequencies, the pure 
amplification dominates with phases distributed around and —2n. The dip observed in the 
IMD3 curves around 1 /t^ therefore corresponds to a transition between these two regimes, 
where the two contributions tend to cancel each other as they are of opposite signs. This is also 
confirmed by the TT-phase shift at the dip frequency that is seen on figure |3(b)| It can also be 
noticed that when the SOA gain gets weaker, as for the solid blue curve (4,a.v = 75 mA), then 
the amplification process always dominates over the CPO effect, even if they are probably quite 
balanced at low frequencies, explaining the low level of IMD3 in this frequency domain. 

5. Spurious-Free Dynamic Range 

The spurious-free dynamic range (SFDR) is a key parameter to evaluate the performance of 
a microwave photonic link in a radar system. It is determined by the noise floor PNoise, and 
the third order intercept point IP3. This point corresponds to the extrapolated fundamental RF 
output power Pi such as Pi = IMD3. The noise floor is defined as the electrical noise power 
contained in a 1 Hz electrical analysis bandwidth. The values of PNoise we used in this paper 
were both measured and theoretically calculated for our SOA as a function of the bias current 
and the RF frequency 1231 . Finally, the SFDR is obtained according to 1(121 : 



and corresponds to the situation where the spurious (here the IMD3) equals the noise floor, as 
illustrated in|4] 
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Fig. 4. SFDR determination for 4/^,, = 250niA, Pin = lOmW, and IMD^" = Ji{m)J2{m). 
In blue: low modulation frequency (<C 1 / tO; in red: high modulation frequency 1 / T^); 
in green: frequency around the IMD3 dip 1/ is)- 



We calculated both the SFDR and the third order intercept point IP3 in the modulation 
frequency range from 40kHz up to 30GHz, and for various injected currents. Once again, 
we compared the situation of a perfectly linear modulation, namely without intermodulation 



products at the input of a SOA, with the case of an actual Mach-Zehnder modulation, with 
intermodulation products given by the Bessel functions expansion. These two situations are 
represented in figures [5(a)| and [5(b)| respectively in dotted and solid lines. 




Fig. 5. (SFDR (a) and 3 order intercept point IP3 (b), as a function of the modulation 
frequency, for various injected currents, and in the two configurations of a perfectly linear 
modulation and of a realistic Mach-Zehnder modulator (respectively in dotted and solid 
lines). 



A first conclusion is that in both cases, over the full spectral range of interest for slow light 
applications (between 0.1 and lOGHz), and for all considered bias currents, the SFDR of the 
simulated link scales from 105 to 125dB/Hz^/^. These values reasonably match radar system 
requirements. Moreover, for modulation frequencies below or in the range of I/Tj, the pres- 
ence of an initial IMD3 at the input of the SOA can even improve the SFDR of the link. This 
phenomenon is widely known among the analog optical transmission community, where the 
nonlinear gain transfer function of a SOA is used to linearize a Mach-Zehnder modulator and 
hence reduce the nonlinear distortion effects 12411251 . However, our theoretical analysis of this 
effect in the framework of microwave photonics, and its interpretation by the input IMD3 being 
in antiphase with the CPO gain gratings (see section|4|, cast new light on this effect, and gives 
it different perspectives when slow light applications are considered. Moreover, it is also impor- 
tant to notice that the frequency domain where a dip is observed in the IMD3, and consequently 
a peak in the SFDR, roughly corresponds to the domain where slow-light effects are the most 
efficient in the SOA, i.e. where the CPO-induced phase shift is maximal. This last point could 
be of importance from the operational point of view. 

6. Conclusion 

We reported in this paper a generalized and predictive model that enables to derive the high 
order harmonic generation and the intermodulation products inside a SOA. This model takes 
into account the variations of the gain saturation parameters along the propagation axis into the 
SOA, and uses a rigorous expression of the gain oscillation harmonics. We also showed that 
provided a basic small signal gain measurement is conducted, this model only relies on adjust- 
ment parameters which are material constants, ensuring the predictive capability of the model. 
This enabled us to derive the SFDR of a SOA-based microwave photonic link for various exper- 
imental conditions. In particular, we investigated the cases where either a perfectly linear or an 
actual modulator is used at the input of the link. We came to the conclusion that the nonlinear 
distortion induced by the input modulator could even be favorable to the spurious-free dynamic 



range of the link, especially in the spectral range where slow light induced phase shifts are the 
most important. Within the field of analog optical transmission, the usual interpretation of this 
linear dynamic range improvement is the compensation of the nonlinearities of the sinusoidal 
transfer function of the Mach-Zehnder modulator by the nonlinearities of the gain transfer func- 
tion of the SOA. However, by considering slow light applications, we proposed a more in-depth 
interpretation based on the CPOs in the semiconductor amplifier and including RF frequency 
dependence, which we believe is therefore more powerful and realistic than the usual one. Fi- 
nally, the results obtained in this work are encouraging since they show that the nonlinearities 
brought by a commercial SOA, when operated in the slow light regime, are already acceptable 
for the less demanding radar applications. 
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